Paramagnetic Spin Correlations in CaFe2As2 Single Crystals 
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Magnetic correlations in the paramagnetic phase of CaFe2As2 (Tat = 172 K) have been examined 
by means of inelastic neutron scattering from 180 K (~ l.OSTiv) up to 300 K (1.8rjv). Despite the 
first-order nature of the magnetic ordering, strong but short-ranged antiferromagnetic (AFM) cor- 
relations are clearly observed. These correlations, which consist of quasi-elastic scattering centered 
at the wavevector Qafm of the low-temperature AFM structure, are observed up to the highest 
measured temperature of 300 K and at high energy transfer (hu > 60 meV). The L dependence 
of the scattering implies rather weak interlayer coupling in the tetragonal c-direction corresponding 
to nearly two-dimensional fluctuations in the (aft) plane. The spin correlation lengths within the 
Fe layer are found to be anisotropic, consistent with underlying fluctuations of the AFM stripe 
structure. Similar to the cobalt doped superconducting BaFe2As2 compounds, these experimen- 
tal features can be adequately reproduced by a scattering model that describes short-ranged and 
anisotropic spin correlations with overdamped dynamics. 



PACS numbers: 74.70.-b,75.30.Et,78.70.Nx 

I. INTRODUCTION 

The appearance of superconductivity (SC) in doped 
v4Fe2As2 materials {A — Ca,Sr,Ba) is linked to the sup- 
pression of antiferromagnetic (AFM) ordering found in 
the parent compoundsii"— Other unconventional super- 
conductors share a similar phase diagram, suggesting 
that the AFM spin fluctuations may be responsible for 
pairing of electrons in the SC state. Certainly, the 
AFM fluctuation spectrum itself is directly influenced 
by superconductivity:^^" The appearance of a gap and 
resonance-like feature in the paramagnetic spectrum of 
these compounds below the superconducting tempera- 
ture Tc also closely resembles^ other unconventional su- 
perconductors and highlights the coupling of the iron 
spins with electronic charge carriers. Recent neutron 
scattering results elegantly show that the magnetic reso- 
nance can infer details about the symmetry of the super- 
conducting gap 4 

In order to understand the nature of the superconduct- 
ing pairing, details of the normal state spin fluctuations 
must be first understood. These spin fluctuations are 
expected to be unusual^"~ since magnetic frustration in 
the tetragonal paramagnetic phase leads to an additional 
nematio^ degree-of-freedom due to the weak net mag- 
netic coupling of nearest-neighbor Fe sublattices. The 
AFM ordering in parent A¥e2^^2 compounds below T/y 
is observed to occur either simultaneously with or after 
a structural transition from tetragonal to orthorhombic 
at Ts. 

The magnetic excitations in AYe2^^2 in stripe 
AFM ordered parent compounds have been extensively 
studiedi^"— and indicate itinerant spin waves whose ex- 
citation spectrum can be adequately described by large 



nearest (Ji) and next-nearest-neighbor (J2) in- plane ex- 
change constants, a substantial interlayer coupling (Jc), 
and strong Landau damping 7. 

The primary role of doping is to suppress both the 
AFM and structural transitions, resulting in strong AFM 
spin fluctuations. In the context of establishing the con- 
nection between spin fluctuations and superconductivity, 
it becomes imperative to examine the evolution of these 
spin excitations across the entire phase diagram of doped 
v4Fe2As2. There have been several theoretical and ex- 
perimental studies on the effects of charge doping on the 
spin excitations in the superconducting compoundsi^^— 
In general, these studies show that the normal state spin 
fluctuations are quasi-two-dimensional (2D) and strongly 
damped while remaining peaked at Qafm, the magnetic 
ordering wavevector of the stripe AFM state. It is un- 
clear if the quasi-2D nature and strong damping are fea- 
tures that appear only with sufficient doping, or whether 
these are common features of the spin fluctuations in the 
tetragonal phase. 

In order to address these questions, we examine the 
spin correlations that occur in the parent compound 
CaFe2As2 above the AFM ordering temperature T/v. 
Inelastic neutron scattering data on CaFe2As2 single- 
crystals show that above the simultaneous first order 
transitions at T/v = Ts = 172 K,^^ the spin gap col- 
lapses and spin wave scattering centered at the stripe 
ordering wavevector Qafm is replaced by short-ranged 
and quasi-elastic AFM correlations that extend up to at 
least 60 meV. Just above Tn (at T = 180 K) the low en- 
ergy magnetic response is quasi-elastic with anisotropic 
in-plane correlations. We find the in-plane correlation 
length to be — 8 A along the orthorhombic a-axis 
and —6 A along b. Weak modulations of the scat- 
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FIG. 1: Summary of scans performed at T = 10 K, 140 K (empty 
symbols) and 180 K (filled symbols) on HB3 {Ef = 41.2 meV) and 
HB-IA with spectrometer configurations described in the text, (a) 
shows (hhL) plane in reciprocal space where the scans at HB3 were 
performed, (b)-(f) shows the various cuts investigated, as indicated 
in (a), (g) and (h) show the scans performed in the (hOL) plane 
using HBIA. No diffuse magnetic scattering was observed in the 
(hOL) plane. The magnetic signal is only observed centered at 
wavevectors Q = Qafm- 



tered intensity are also observed along the c-axis indi- 
cating a two-dimensional character to the paramagnetic 
fluctuations. In general, spin correlations weaken and 
broaden further in momentum and energy with increas- 
ing temperature, but are still observed up to the highest 
measured temperature of 300 K. These observations can 
be explained in the context of spin dynamics overdamped 
by particle-hole excitations. In particular, we use a phe- 
nomenological theoretical model with in-plane and inter- 
plane magnetic anisotropy to consistently fit our data for 
all temperatures, obtaining the ratios J1/J2 — 0.55 and 
Jc/J2 — 0.1. We find that the spin fluctuations in the 
paramagnetic phase of the parent compound bear a close 
resemblance to the paramagnetic fluctuations in the su- 
perconducting compositions. 

This article is laid out as follows. In section II below, 
the experimental conditions under which the experiment 
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FIG. 2: Temperature evolution of the neutron scattering signal 
measured on HB3 at huj = 10 meV. (a) Background estimate 
measured away from Qafm a-t (002) (solid symbols) and (003) 
(open symbols) showing no anomaly at T^y. Solid lines are lin- 
ear fits to the temperature dependent intensity, (b) Intensity at 
Qafm = (1/2, 1/2, 1) and (1/2, 1/2, 0). The solid line is the 
non-magnetic background estimate obtained from averaging the 
fits at (002) and (003), shown in panel (a). 



were performed and the sample details are presented. 
The data analysis and results are presented in section 
III. Finally, a discussion and a summary are given in sec- 
tion IV. 



II. EXPERIMENTAL PROCEDURES 

Inelastic neutron scattering measurements were per- 
formed on a single crystal mosaic (^400 small single- 
crystal samples) of CaFe2As2 with a total mass of 
^ 2 grams that are co-aligned to within 1.5 degrees 
fuU-width-at- half- maximum (FWHM). The preparation 
methods of the single-crystals have been described 
elsewhere.^* Data were collected using the HB3 and 
HBlA triple-axis spectrometers at the High Flux Iso- 
tope Reactor at Oak Ridge National Laboratory and 
the MAPS chopper spectrometer at the ISIS facility at 
Rutherford Appleton Laboratory. HBB was operated in 
relaxed resolution for measurement of the diffuse scat- 
tering signals in the paramagnetic phase, with flxed final 
energy (£'/) configurations, Ef = 14.7 meV and 41.2 
meV, and 48'-60'-80'-120' coUimation. The sample was 
mounted in a closed-cycle refrigerator and oriented for 
scattering in the tetragonal {hhL) plane. HBIA was op- 
erated with fixed incident neutron energy of 14.7 meV 
and 48'-40'-40'-136' coUimation and the sample mounted 
in the (hOL) plane. The MAPS experiment was per- 
formed at T =180 K, with an incident energy of 100 meV 
using the same sample aligned with the c-axis along the 
incident beam direction. 

To avoid confusion, the data is exclusively presented in 
tetragonal units and we define Q — ^{hi + kj) + ^Lk 
as the momentum transfer indexed according to the 
lA/mmm tetragonal cell with lattice parameters a = 3.88 
A and c = 11.74 A at 300 K. The vectors i, j and k are the 
fundamental translation unit vectors in real space. For 
comparison with the AFM low temperature orthorhom- 
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FIG. 3: (Color online). L and h-dependence of the scatter- 
ing at huj = 10 meV measured on the HB3 instrument with 
Ef =41.2 meV. The scans are performed along the (1/2, 1/2, L) 
and (h, h, 3) directions for temperatures T =300 K [(a) and (b)] 
and T =180 K [(c) and (d)]. The solid lines in (a) and (c) corre- 
spond to fits to the dynamical susceptibility described in Eq. (|5]l. 
The solid lines in (b) and (d) are guide to the eye, and based on 
Lorentzian fits to the data. The dashed line is an estimate of back- 
ground scattering. Panels (c) and (d) show the sharper magnetic 
scattering at 180 K. 



bic (o) structure, we note the following relations between 
the Miller indices of the two phases, h^{Ho + Ko)/2, 
k={Ho — Ko)/2, and L=Lo. For convenience, we some- 
times use the reduced momentum transfer q=Q — Qafm 
in reciprocal lattice units (rlu) where Q is the momen- 
tum transfer to the sample and Qafm=(^Oi ^o, ^o) the 
reciprocal lattice vector which defines the AFM low tem- 
perature zone center. Typical AFM wave vectors studied 
are (1/2, 1/2, L) with L = odd. 

After suitable subtractions of the non-magnetic (back- 
ground) scattering, the observed magnetic inelastic neu- 
tron scattering data is cast in terms of the dynamical 
structure factor S{CI,lu), which is related to the imagi- 
nary part of the dynamic spin susceptibility x"(Q,a;) via 
the fluctuation-dissipation theorem. 



5(Q,w) = 2(ro 



47r^| 1 - e-^"/*^^ 



(1) 



where (rp) — 290.6 mbarns Sr~^ is a conversion fac- 
tor to bring the intensity into absolute units of mbarns 
meV~^Sr~^ f.u."-'^ (Sr=Steradian, f.u.=Formula Unit) 
and F(Q) is the magnetic form factor for the Fe^+ ion. 



III. ANALYSIS AND RESULTS 

A. Survey of reciprocal space 

Fig. [T]shows the (hhL) plane in reciprocal space where 
the low energy measurements were performed. It also 
shows several Q-scans taken at temperatures below (10 
K and 140 K) and above Tjv (180 K) the concomitant 
structural and Neel ordering temperature T/v = = 172 
K. Below T/v, magnetic Bragg peaks appear at Qafm = 
(1/2, 1/2, L) positions with L = odd that describe the 
ordered AFM stripe structure. Figs, mb)-(f) show var- 
ious cuts, as indicated in Fig. [Ua), through the (hhL) 
scattering plane at a finite energy transfer of 10 meV and 
Ef — 41.2 meV. Similar scans performed with Ef = 14.7 
meV show qualitatively the same results. Below T/v, 
sharp excitations are observed at Qafm that originate 
from very steep spin wave excitations in the ordered 
stateJ^ii^ The difference in sharpness of the spin wave 
peaks in the [h, h, 0] and [0, 0, L] directions is due to 
the anisotropy in the spin wave velocity, as discussed 
in Ref . [li| . When the sample is warmed up above T/v, 
strong intensity remains at Qafm position with much 
broader lineshapes indicating short-ranged AFM spin 
correlations. The scans shown in Fig. [1] are strongly 
influenced at higher angles by the presence of aluminum 
phonon scattering from the sample holder and low angle 
scattering from the direct beam. Both contributions lead 
to very high background levels and limit the range of the 
Q-scans. 

Since the stripe AFM ordering may be frustrated in the 
tetragonal structure, we searched for evidence of spin cor- 
relations at other wavevectors in addition to the strong 
components of the diffuse excitations near Qafm- No 
magnetic diffuse scattering was observed along various 
symmetry directions in the (hOL) plane. In particular, no 
evidence of magnetic scattering was seen at the wavevec- 
tor (1,0,L = even) corresponding to Neel (C-type) AFM 
fluctuations (see Figs. [TJg) and (h)). In the (hhL) plane, 
there are indications of weak peaks in the extended Q- 
scans at wavevectors other than Qafm which may arise 
from additional magnetic modulations in the paramag- 
netic phase. For example, very weak peaks can be ob- 
served at (002) and (003) (see Figs, [^b), (d), and (e)) 
which would correspond to the presence of ferromagnetic 
correlations and A-type magnetic correlations (ferromag- 
netic within the layer, AFM between layers), respectively. 

In order to understand the development of correlations 
at Qafm and address the potential existence of addi- 
tional magnetic modulations (002) and (003), the tem- 
perature dependence was measured at various points in 
the {hhL) plane at an energy transfer of 10 meV. Fig. 
Ufa) shows the temperature evolution of the scattered 
intensities at (002) and (003). The intensities show no 
anomaly at T/v, but rather the intensity increases linearly 
with temperature as expected for a phonon background. 
This background is observed throughout the scattering 
plane and partly arises from aluminum phonon scatter- 
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FIG. 4: (Color online) The dynamic magnetic susceptibility of 
CaFe2As2 as a function of energy at Qafm = (1/2, 1/2, 3) for 
T = 140 K (blue solid squares), T = 180 K (grey solid circles) and 
T = 220 K (red solid diamonds). Above the magneto-structural 
transition Tjv = Tg = 172 K, a broad magnetic spectrum is ob- 
served as quasi-elastic response near Qafm- Data taken at T = 180 
K were fit to a Lorentzian form given in Eq. ((2]| convoluted with 
the instrumental resolution (black solid line) . The Lorentzian half- 
width Ft at T = 180 K is 10 meV. At T = 220 K, we estimate the 
energy linewidth to be ~ 13 mcV using the expression Tt = 7(^)'^ 
and the fitted value of 7 (temperature independent Landau damp- 
ing defined in the text) and that of the correlation length at 220 
K. The calculated Lorentzian scattering at T = 220 K is shown as 
a red solid line. In contrast, sharp spin waves having an energy gap 
A of ~7 meV are observed in the ordered phase at T = 140 K. 



ing from the sample holder, with little or no magnetic 
contribution. We can use the intensity at (002) and (003) 
as an estimate for this phonon background. In Fig. EJb), 
the temperature dependence at Qafm — (1/2, 1/2, 1) 
and also (1/2, 1/2, 0) is shown. The intensity at Qafm 
increases from low temperatures as expected for the in- 
creasing Bose population factor of the low lying spin wave 
modes. At T^, there is a sharp drop in the intensity con- 
sistent with the first-order transition to the paramagnetic 
state. The paramagnetic intensity decreases above T^r 
and is nearly at background by 300 K. At (1/2, 1/2, 0) 
the intensity is nearly at background level below Tjv, in- 
creases sharply at the transition, and decreases slowly at 
higher temperatures. 

The L— dependence of the paramagnetic scattering is 
consistent with weak antiferromagnetic correlations be- 
tween layers atop a constant magnetic background, as 
illustrated in Fig. [2][b). Figs. [31(a) and (c) show the L- 
dependence of the scattering &i fiijj — IQ meV and along 
the (1/2, 1/2, L) direction for T = 300 K and 180 K, 
respectively. The L-dependence displays sinusoidal vari- 
ation with maxima at odd values of L. The lines shown 
correspond to fits to the dynamic susceptibility and will 
be described in detail below. The intensity drops by a 
factor of 2 between 180 K and 300 K indicating the grad- 
ual evolution of the system to less correlated quasi-2D 
spin fluctuations similar to the reduction in intensity for 
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FIG. 5: (Color online). Magnetic excitations in CaFe2As2 mea- 
sured on the MAPS spectrometer with an incident energy of 
Ei = 100 meV at T = 10 K (left panel) and 180 K (right panel). 
The data shows magnetic intensity as a function of the [h, h\ di- 
rection and the energy transfer after averaging over the transverse 
[h, —h\ direction in the range 0.4< h <0.6. Given the fixed crystal 
orientation with incident beam along L, the L component of the 
wave vector varies with the energy transfer as indicated. Excita- 
tions below Tjv are consistent with steep spin waves and diffuse 
magnetic excitations are observed above Tjv at T = 180 K. 



cuts along the [/i, /i, 3] direction shown in Figs. [3] (b) 
and (d). The temperature dependence of Al phonons is 
primarily responsible for the increase in background be- 
tween 180 K and 300 K. However, the constant magnetic 
background itself is also weakly temperature dependent 
as inferred from Fig. [2] (a). We return to this in section 

Unci 

In all, the surveys of magnetic scattering intensities 
above T^v in the (hhL) and (hOL) planes indicate that 
the AFM spin correlations are restricted to the vicinity of 
Qafm, the wavevector of the stripe ordered phase. One 
essential difference in the paramagnetic phase is an in- 
crease in the c-axis anisotropy and tendency towards 2D 
spin fluctuations, as indicated by the weakly modulated 
rod of scattering along L. This is entirely analogous to 
the behavior of AFM spin fluctuations in the optimally 
doped superconductors, where interlayer correlations are 
very weak^i^. 



B. Energy dependence 

Fig. m depicts the dynamical structure factor 
x"(QAFM,w)/?iw at Qafm = (1/2, 1/2, 3) and huj < 22 
meV for temperatures above and below Tpf. In order to 
obtain x" from the raw data, a non-magnetic background 
was estimated by measurements at Q = (0.35, 0.35, 3) 
and subtracted. Below T/v, the low energy magnetic 
spectrum consists of spin waves with a sizeable spin gap 
of 7meV (see Ref.[il). Just above T^r, the gap collapses 
and the spin wave scattering is replaced by gapless, dif- 
fusive excitations. At 180 K, the diffusive excitations can 
be fit to a quasi-elastic Lorentzian form. 



x"(Qafm,^) 
huj 
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with an energy linewidth of T-r = 10 ± 1 meV. The 
parameter A is an arbitrary intensity scale factor. As 
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FIG. 6: (Color online). Constant energy slices (Ai? = ±5 meV) 
through the excitation spectrum of CaFe2As2 in the {h, k) plane 
at T = 10 K and 180 K, as observed on MAPS. Energy slices 
are chosen to correspond to odd values of L. Intensity shown is 
in absolute units (mbarn Sr~^ meV~^ per formula unit). Below 
Tjv, well defined spin waves are observed around Qafm (see also 
Ref . [T^ ) . Above Tm , strong but short-range magnetic correlations 
remain around Qafm a-nd extend up to at least 60 meV. Solid lines 
in panel (c) show the directions along which the cuts in Fig. [7] are 
taken, with (+) designating longitudinal cuts and (-) transverse 
cuts. 



the temperature is increased, the Lorentzian half- width 
grows rapidly. Kt T — 220 K, the spectrum weakens 
considerably with temperature and appears flat within 
the energy range measured, thus the Lorentzian widths 
become large and poorly defined. This rapid increase in 
the quasi-elastic linewidth with temperature is explained 
below. 



C. Spectrum of paramagnetic spin fluctuations 
near Qafm 

At temperatures just above T/v (T = 180 K, T/Tn = 
1.05) we used the MAPS spectrometer to perform a de- 
tailed survey of the spin fluctuations in the paramagnetic 
phase in the vicinity of Qafm- For detailed modeling, the 
MAPS measurements were normalized in absolute scat- 
tering units of mb Sr"""^ meV~^ f.u.~^ by comparison to 
a vanadium standard. Data were collected at 180 K and 
also in the AFM ordered phase at 10 K (with an inci- 
dent energy of 100 meV). As in previous work,^® we use 
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FIG. 7: Longitudinal (-|-) and transverse (-) constant-energy cuts 
measured around Q = Qafm=(1/2, 1/2, L) on the MAPS spec- 
trometer. The cuts arc taken at hu> =12, 39 and 60 meV, which 
correspond to L = 1, 3 and 5, respectively. Solid lines are best fits 
of Eq. ^ to the data. Dotted lines indicate the fitted background, 
and the dashed lines show the instrumental resolution function. 



the MSLICE program^^ to visualize the data and to take 
one and two dimensional cuts through main crystallo- 
graphic symmetry directions for subsequent data analysis 
with the TOBYFIT suite of analysis programs described 
below. ■^^ Where possible, symmetry equivalent cuts and 
slices were added to improve statistics. Fig. [5] shows 
slices of the neutron intensity along the [/i, h] direction 
as a function of energy transfer after averaging over the 
[/i, — h] direction. Below Tjv, the data in the left panel of 
Fig. [5] show a steep plume of intensity centered at Qafm 
arising from AFM spin waves. The sizeable exchange 
coupling along c leads to variations in the structure factor 
along L, observed as energy-dependent intensity oscilla- 
tions that are peaked at the AFM zone centers; hijj — 12 
meV {L=l), 39 meV (L=3), and 60 meV (L=5). Analy- 
sis of the AFM spin wave spectra is described in detail in 
Ref. [ll]. Above Tjv, the right panel of Fig. [5] indicates 
that the magnetic spectrum is much broader in Q, and 
energy dependent oscillations are much less pronounced, 
confirming short-ranged spin correlations within the Fe 
layer and a weakening of interlayer correlations. 

Fig. ini shows the neutron intensity for several constant 
energy slices at T = 10 K and 180 K. The energies are 
chosen such that L = 1, 3, and 5, in order to measure the 
spin fluctuations in {h, fc)-planes containing Qafm (i-e. 
measuring the spin correlations within the Fe layers) . Be- 
low T/v, Figs. [n](a)-(c) show three constant energy slices 
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FIG. 8: Net magnetic neutron scattering intensity as a function of 
energy measured on MAPS around Qafm = (1/2, 1/2, L) after 
subtraction of a non-magnetic background at Q = (0.15, 0.5, L). 
The integration ranges used in MSLICE are such that Ah = ±0.05. 
Lines are the best fits of Eq. js} to the data. 



through the AFM spin wave cone centered Qafm- Above 
Tjv, Figs. |6](d)-(f) again demonstrate that the sharp spin 
waves are replaced by very broad scattering in the para- 
magnetic phase. The paramagnetic fluctuations persist 
up to at least 60 meV and are notably anisotropic with 
an elliptical shape. This anisotropic scattering is charac- 
terized by transverse [h,—h] cuts through Qafm being 
broader than the longitudinal [/i, h] cuts, as illustrated 
in Fig. [7] where the solid lines are the resolution convo- 
luted fits of Eq. ([5]) to the experimental data, described 
below. This anisotropy implies that the paramagnetic 
fluctuations are described by two different spin correla- 
tion lengths within the Fe layer. An in-plane anisotropy 
of this type is allowed in the I4/mmm tetragonal cell, i.e. 
4-fold symmetry has not been broken. The origin of this 
anisotropy is discussed below. 

Fig. m shows the energy dependence of the MAPS neu- 
tron intensity at Qafm = (1/2, 1/2, L) at T = 180 
K after subtraction of background estimated at Q = 
(0.15, 0.5, L). In the geometry of our time-of-fiight 
experiment, L varies with energy transfer, however, L- 
dependent oscillations here and in Fig. [5] are only weakly 
present. After correcting for the dependence of the mag- 
netic form factor, the data was fit to the full susceptibility 
function described in Eq. [5] after convolution with the in- 
strumental resolution function. The model describes the 
MAPS data quite well, including weak L-dependent os- 
cillations, and gives relaxational half-width Tt = 10±1 
meV, which is consistent with the IIB3 measurements 
shown in Fig. 

Fig. ini shows IIB3 scans along the [h, h] and L- 
directions through Qafm = (1/2, 1/2, 3) and ?ia; = 5 
meV at several temperatures. Similar to Fig. [1] the 
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meV. Solid lines are the best fit to the data. Both the h- and L 
scans were fit to Eq. JSjl (as shown in Fig. |3]over a wider range of 
h- and L). The data at T = 140 K was fit to a spin wave model 
(as described in Ref. [Tsll ). and is shown here for comparison. 



cuts below Tn show sharp spin wave scattering. The 
measurements at several temperatures above Tpf indicate 
a gradual weakening of the scattering at Qafm as the 
temperature is raised, reminiscent of Figs. [3] and |4l As 
noted above, the magnetic scattering along the [h, h] di- 
rection through Qafm is much sharper, as shown in Fig. 
[3] (d), indicating much longer ranged spin correlations 
within the Fe layer than between layers. We see now that 
the weakening originates from temperature-dependent 
broadening of the scattering in reciprocal space, indicat- 
ing a gradual reduction in the spin correlation length 
with temperature. Similar to Fig. [31 the L— dependence 
of scans at Qafm =(1/2, 1/2, 3) and hut — 5 meV were fit 
to the susceptibility in Eq. ([5]) for several temperatures 
up to 300 K and are shown as solid lines in Fig. [HUb). 
The description of this fitting and the interpretation of 
the results is now described. 



D. Theoretical Model 

In this section, we describe the theoretical model used 
to fit the experimental data. At temperatures below 
Ttv, the imaginary part of the generalized susceptibil- 
ity x"(Q,w) for CaFe2As2 is well described by damped 
spin wave modelsj^^"— Above Tn, our neutron scatter- 
ing data reveal several features that have allowed us to 
develop a fairly detailed model of the paramagnetic exci- 
tations. First, the excitations are diffusive in character, 
at least for the measured energy scales (< 60 meV). Sec- 
ond, the excitations are nearly two-dimensional, with a 
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weak modulation of the scattered intensities along L. Fi- 
nally, broad scattering is observed within the Fe layers 
that is anisotropically distributed around Qafm- We are 
thus forced to define two in-plane spin correlation lengths 
in order to fully describe the quasi-2D spin fiuctuations. 
Considering these facts, we start by writing the to- 



tal dynamical magnetic susceptibility per formula unit 
as that of a two-dimensional metallic AFM in the para- 
magnetic phase,— >21 which we have extended to include 
the effects of an in-plane anisotropy of the magnetic in- 
teractions. 



J 



AFM 



qii, 



{ql + mxqy)a^ + {£,T/a) ^ - ihu/'j 



r 



(3) 



The parameter defines the magnetic correlation length 
at temperature T in units of A, xo is the staggered sus- 
ceptibility, and 7 denotes the damping coefficient origi- 
nating from the spin decay into particle-hole excitations. 
The dimensionless parameter rj is used to represent the 
anisotropy of the in-plane correlation lengths. The vector 
Qll = {Qx, Qy) = ^{h, k) is a 2D vector defined relative to 
Q2i>j^ = "^{1/2, 1/2) and a = y/2ao is the lattice con- 
stant of the tetragonal cell with oq the nearest neighbor 
Fe-Fe distance. We can define two correlation lengths 
in directions longitudinal (qx = % denoted (-1-) scan) 
to Qafm '^i^h a value = fT(l + f )^/^ and trans- 
verse {qx = —qy denoted (-) scan) to Qafm '^^^^ value 
= <^t(1 - f)^/^- Thus, we must have < < 2 

I 



in order for these correlation lengths to be meaningful, 
which we found indeed to be the case. 

The effect of weak interlayer AFM coupling Jc on total 
dynamic susceptibility can be written as 



X ^QAFM + q,'^) = X2D(QAFAf + q|h'^) + 2>/cSin^(g2c/4) 

(4) 

where the scattering vector is defined as Q = Qafm + 
q|| -|- g^k. This expression can be motivated using the 
1/iV approach for the action of the collective magnetic 
degrees of freedom employed in Refsi^i^. Taking the 
imaginary part of Eq. we obtain the expression 



x"(Q,^) 



{<f + V<lxqy)a^ + (^)-2 + + cos(7ri)) 
I 



2 ' 



(5) 



TABLE I: Best fit parameters of Eq. Q obtained self consis- 
tently from both HB-3 and MAPS at temperature T =180 K. The 
parameter xo was solely obtained from the MAPS data which were 
normalized in absolute units. 
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Xo 


V 




(meV) 


(A) 


(Ats^ meV"' f.u."') 






43± 5 


7.9±0.10 


0.2±0.05 


0.55±0.36 


0.20±0.02 



which we used to fit our inelastic neutron scattering data. 
Here, the parameter rjc = JcXo controls the strength of 
the interlayer spin correlations and we have substituted 
q^ = 2'k{L — l)/c. A physical motivation for this ex- 
pression and an interpretation of these parameters in the 
context of the iron arsenides is given later in section IV. 

As we are only able to define -q from a single MAPS 
data set at 180 K, we make the assumption that the pa- 
rameters 7j Jc and rj in Eq. ([S]) do not depend on tem- 



perature, making the only temperature dependent pa- 
rameter in the model. However, one can introduce effec- 
tive temperature dependent parameters Xt = Xo(^)^ 
and Ft = 7(^)^- With these definitions, the dynam- 
ical susceptibility along the longitudinal direction at, 
q — (y^-^, L — Ij, or in the 2D limit where Jc = 
can be written as. 



X2DiQ'AFM 



q,w) = 



r|(i 



(6) 



At g = 0, this model reduces to the quasi-elastic 
Lorentzian lineshape used to describe the susceptibility 
inEq. ©. 

In the context of the iron arsenides, Eq. (|6]) was first 
used successfully by Inosov et ali^ to describe the Co- 
doped BaFe2As2 superconducting compound. In their 
work, Inosov et al^ assumed a mean-field behavior for 



the correlation length (T — T, 



>l/2 



Here, given the 
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TABLE II: Temperature evolution of the correlation length 
(in A) for several incident energies, as measured on HB-3. 
For a given temperature, appears largely independent of 
energy, as expected. The final column, containing the corre- 
lation length averaged over all energies, is the value used in 
subsequent model calculations. 





5 mcV 


10 meV 


15 meV 


< Ct > (A) 


180 K 


8.1±0.2 


7.9±0.1 


7.8±0.3 


7.9±0.1 


220 K 


7.0±0.2 


7.2±Q.3 


6.9±0.2 


7.0±0.1 


250 K 


6.0±0.1 


6.0±0.2 


6.2±0.2 


6.1±0.1 


300 K 


4.9±0.1 


4.9±Q.2 


4.9±0.2 


4.9±0.1 



strong first-order character of the magnetic transition in 
the parent compound CaFe2As2, we do not assume any 
particular form for temperature dependence of S,t ■ 

For detailed comparison of model susceptibility in Eq. 
([5|) to the data, we use standard inelastic neutron scat- 
tering computer programs that account for resolution 
effects. Specifically, the high energy data collected on 
MAPS was analyzed using TOBYFIT^^ while the low 
energy response measured on HB3 (up to 15 nieV) was 
analyzed using RESLIB.^^ 

Our protocol was to establish the value of the temper- 
ature independent parameters using data sets from both 
HB3 and MAPS at 180 K where scattering features are 
more sharply defined. Data fits at subsequent tempera- 
tures would then only require the refinement of a single 
parameter, ^t- Due to the interrelated nature of the fit- 
ting parameters, self-consistency checks were performed 
for model fittings on the HB3 and MAPS data sets. The 
in-plane anisotropy parameter rj was obtained indepen- 
dently from the MAPS data by fitting the longitudinal 
{+) and transverse (-) cuts, as depicted in Fig. [S] The ob- 
served in-planc anisotropy is then determined from these 

fitted widths according to ?] = ^J^. The temper- 

ature independent damping parameter, 7, was obtained 
from both the HB3 and MAPS energy scans at T = 180 
K where we have the best signal to noise ratio. The cor- 
responding fits are illustrated in Figs. 0] and [71 The fits 
were performed by constraining the value of 7 = Tt{^)'^ 
such that the Lorentzian width Tt remains fixed at 10 
meV (fitting only for ^t). Once 7 was determined at 
T — 180 K, it was held fixed in all subsequent fits. 

To determine rjc at 180 K, we have performed a chi- 
squared analysis of fits in which only rj^ {0 < rjc < 0.5) 
was allowed to vary, and other parameters held fixed. 
The best fits were obtained for 0.9 < rjc < 0.11, as indi- 
cated in Fig. ([3]), which also corresponds to a minimum 

value. Outside this range, drastically increases and 
the fits are altered. The parameter 77c in Eq. ([5]) which 
describes the interlayer coupling was found to be ~ 10% 
at T = 180 K. 

The parameters xo and were found to be strongly 
correlated. The prescription was used to fix (effec- 



tively fixing the width) and then vary the scale factor xo 
to get the best fit. This procedure was performed for con- 
stant energy cuts at several energy transfers on both HB3 
and MAPS data sets at T = 180 K. Upon convergence 
to reasonable values of xo and ^t, the whole procedure 
was repeated for consistency. Once the best fitted value 
of Xo was obtained at 180 K, we kept it fixed at that 
same value for all other temperatures and fit the data by 
varying only ^t- 

The MAPS data was normalized to vanadium stan- 
dard, allowing us to report xo in absolute units. In con- 
trast, the HB3 data was not in absolute units. Thus, 
the corresponding best fitted scale factor from HB3 gives 
a value that is proportional to the true xo- Therefore 
the Xo value reported here are obtained solely from the 
MAPS data. In all, we arrived at the following final 
set of parameters at T = 180 K; 7 = 43 ± 5 meV, 
Xo = 0.20 ± 0.05 Ai| meV-if.u.-i, Ct+ = 8 ± 1 A and 
= 6± 1.5 A. Fits are shown as solid lines in Figs. |31 
m [7] and [5] From these values, we estimate the average 
in-plane anisotropy parameter (ry) ~ 0.55±0.36, and thus 
kept 77 fixed at this value in subsequent model fits to the 
higher temperature data. The corresponding value of the 
static correlation length obtained from MAPS is 7±2 
A, in good agreement with the low energy HB3 value of 
7.9±0.1 A. A summary of all of the model parameters at 
T = 180 K are listed in Table HI 

Fig. [nja) shows the temperature dependence of the 
HB3 constant energy cuts along the [1, 1,0] direction at 
5 meV for several temperatures up to 300 K. These data 
were fit to determine the values of at different tem- 
peratures. The correlation length decreases by a factor of 
~2 from 180 K to 300 K. Similar fits were performed at 
constant energy cuts of 10 meV and 15 meV (not shown), 
giving comparable results. The results are summarized 
in Table HH 

We note that the set of parameters in Tables U and HI] 
describe the full data set at all temperatures and energy 
transfers exceptionally well, as can be seen by reviewing 
all of the fitted curves in Figs. |3 H H [5] and HI For 
example, the line through the data in Fig. |3] at T = 220 
K is not a fit, rather it is a model calculation based on 
the fixed parameters in Tables [T] and [III This universal 
agreement with the data at all measured energies and 
temperatures is a strong endorsement for the validity of 
the nearly AFM model. 

IV. DISCUSSION AND CONCLUSIONS 

A. Discussion 

Despite the strong first order transition and the pos- 
sible influence of magnetic frustration in the tetragonal 
phase, we find that the spin fluctuations observed above 
the AFM ordering temperature in CaFe2As2 are peaked 
at Qafm, the wavevector of the stripe AFM structure. 
We found no evidence of substantial magnetic fluctua- 
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TABLE III: The derived staggered susceptibility Xt, relax- 
ational linewidth Ft, and bulk susceptibility Xbuik, as a func- 
tion of temperature from 180 < T < 300 K. 



Temp. 


180 K 


220 K 


250 K 


300 K 


Xt (m| meV-i f.u "i) 


0.84 


0.66 


0.49 


0.32 


Ft (meV) 


10 


13 


17 


26 


Xbulk (10""' emu mol-i) 


2.41 


2.40 


2.39 


2.37 



tions at any of the other wavevectors studied, indicat- 
ing that the relevant spin fluctuations throughout the 
phase diagram are short-ranged fluctuations of the low- 
temperature stripe AFM ordering. In addition, we find 
weaker interlayer coupling immediately above T/v in the 
parent compound, suggesting that the quasi-2D spin fluc- 
tuations are a property of the tetragonal paramagnetic 
phase in general, and does not necessarily occur only af- 
ter doping. 

The substantial spin fluctuations that persist above Tn 
should have an effect on the bulk transport and magnetic 
properties of CaFe2As2. For example, such spin fluctua- 
tions are a likely source of the high resistivity found for 
CaFe2As2 between and 300 K at ambient pressurei^ 
Large spin disorder scattering at ambient pressure is sup- 
ported by dramatic decrease of the resistivit y^^i'^° and 
disappearance of spin fluctuations^ upon transition to 
the collapsed tetragonal phase under applied pressure. 

The bulk magnetic susceptibility of CaFe2As2 also dis- 
plays a weak, non Curie- Weiss temperature dependence 
above Even though Eq. ([SJ is strictly valid only 

as an expansion for momenta close the AFM ordering 
vector, we can use this expression for the dynamic spin 
susceptibility to estimate the Q = bulk susceptibility as 
Xbulk = x(Q = 0,0. = 0) = Xo/[2^'(l + W2) + i^)-^ + 

For large correlation lengths, the bulk susceptibility 
takes on a temperature independent value of Xbuik ~ 
Xo/[l + 27r2(l + 77/2) + 277c] = 2.43x10-4 emu mor^ 
f.u.^^. This estimate has the same order of magnitude 
as the bulk susceptibility determined from magnetiza- 
tion measurements, ^4 further justifying the validity of the 
model. The shorter correlation lengths found well above 
Ttv will decrease the bulk susceptibility as shown in Ta- 
ble mil although this effect is small for the temperatures 
studied. 

The extrapolation of Eq. ([5]) to other momentum val- 
ues also allows us to estimate the size of the fluctuating 
moment per iron atom 

13h Jx"{Q,u;){l-e-^-/''T)-'dQdLo 
<"^^^2V JdQ • 

The Q-integration is performed in the orthorhombic 
zone (containing one ellipse of scattering) defined by the 
ranges < Q^((3j^) < ^ and < < ^. The factor 
of 1/2 in eq. [7]converts the results from squared-moment 



per formula unit to squared moment per iron. If we 
choose a high energy cut-off of ~200 meV corresponding 
to the upper limit of observed spin wave excitations,— 
we obtain y/JrnF) ~ 0.7 /is per iron which is very close 
to the size of the observed ordered moment 

The spin fluctutations observed above the AFM order- 
ing temperature in CaFe2As2 bear a close resemblance to 
that of doped superconducting BaFe2As2. Based on the 
model fitting, the similarity of the spin fluctuations in the 
parent and superconducting compositions can be com- 
pared more quantitatively. The paper by Inosov et aL~ 
describes the parameters for short-range spin correlations 
in the optimally doped superconducting composition us- 
ing an equation similar to ([6]) . We can estimate the values 
of the model parameters at 180 K for this SC composi- 
tion using mean-field temperature scaling as described by 
Inosov et al. Also, since the model parameters described 
in their work assume isotropic in-plane spin correlations, 
we extended their model by including the anisotropy pa- 
rameter, 77, obtained from our results (i.e. assuming 77 
to be both temperature and composition independent). 
With these assumptions, we estimated the values of the 
model parameters for the superconducting sample with 
the following results at 180 K; the static susceptibility 
at the antiferromagnetic wavevector is Xt = 0.18 
meV"^ f.u.-^, the longitudinal correlation length is £,t^ 
= 11 A, and the Landau damping 7 = 180 meV. This 
can be compared to the values in CaFe2As2 at the same 
temperature as found here; Xt = 0.9 meV~^ f.u.^^, 
^T+ = 8 A, and 7 = 43 meV. The results show that the 
Landau damping parameter is much larger in the super- 
conducting composition, whereas the correlation lengths 
at 180K are similar. In comparing these numbers, it is 
important to also recall the differences between the two 
materials investigated. For instance, in the sample stud- 
ied by Inosov et al.,^^ the alkaline earth metal is Ba and 
no magnetic phase transition is present. On the other 
hand, in the sample analyzed here, the alkaline earth 
metal is Ca and, at 180K, the system is very close to a 
strong first-order magnetic transition. 

We now turn to a discussion of the physical interpre- 
tation of the data. The new feature observed in our 
results is the anisotropy of the in-plane spin fluctua- 
tions in the paramagnetic phase, represented by the pa- 
rameter 77. As discussed above, Eq. ([5]) is based on 
a metallic AFM model and a proper starting point to 
understand this anisotropy is an electronic band model 
for the magnetic susceptibility. The unique band struc- 
ture of the iron arsenides, consisting of circular hole 
pockets at Q=(0, 0) and elliptical electron pockets at 
Qafm=(V2- 1/2), combined with electron-electron in- 
teractions will lead to an anisotropic spin response, 
even in the tetragonal phase. "33'" ^ More physical insight 
can be gained by considering a phenomenological model 
for the collective magnetic degrees of freedom of the 
pnictidesi^"— In these materials, the magnetic lattice 
in each iron arsenide plane can be subdivided in two in- 
terpenetrating antiferromagnetic sublattices, each con- 



taining one iron per unit cell. These two sublattices 
are only loosely coupled, because the local magnetic field 
produced by one sublattice on the other vanishes in the 
completely ordered state. To dominant order, the only 
couplings between the two sublattices that are allowed 
by symmetry are the biquadratic term — g(mi • 1112)^ and 
the momentum-dependent quadratic term riq^qyini ■ 10.2 . 
Evaluation of the spin-spin correlation function using this 
two-sublattice model leads to Eq. ([5]). Within this de- 
scription, the anisotropy of the correlation length results 
from an inter-sublattice coupling. 

This phenomenological model can be motivated by 
either an itinerant^! or a localized moment picture^. 
Indeed, the short correlation lengths encountered in 
the paramagnetic phase suggest a minimum description 
where only the short-range interactions are kept. One 
then has the so-called frustrated J1-J2 model, with an- 
tiferromagnetic nearest-neighbor exchange interaction Ji 
and next-nearest neighbor exchange J2. The anisotropy 
parameter 77 becomes simply the ratio between the two 
exchange constants, 77 = J1/J2, and J2 is related to the 
bare susceptibility xo of Eq. ([5]) through xo = 2/ J2. This 
also allows us to equate the interlayer coupling parameter 
rjc with the ratio 2Jc/J2- 

The exchange ratios in the paramagnetic phase can 
be compared directly to the values in the orthorhom- 
bic AFM ordered phase as obtained from measurements 
of spin waveaiSi^ (see also Fig. |6]) From the observed 
anisotropy of the in-plane correlation lengths, we can de- 
termine J\l J2 = 0.55 ± 0.36 in the tetragonal phase, im- 
plying that Ji is indeed antiferromagnetic. In the AFM 
ordered orthorhombic phase, there are distinct nearest- 
neighbor exchange interactions, J\a and J\h and the ex- 
perimentally determined ratioi^iii (Jia + Jifc)/2J2 ~ 
0.6—1 is similar in magnitude to the exchange ratio 
in the paramagnetic phase and within the error bars. 
Thus, the relative strength of the average exchange in- 
teractions in the Fe plane do not appear to be strongly 
affected by the structural transition. The value of this ra- 
tio Ji/ J2 < 2 places the paramagnetic phase of CaFe2As2 
in a regime of frustrated magnetism within the J\ — J2 
model and consistent with the presence of spin nematic 
correlationsi^"— Taking the fitted value of xo and con- 
sidering the Fe moment to be ~ 1/^b, one can use the 
relation xo ~ 2/J2 to estimate the the effective next- 
nearest neighbor exchange J2 ~ 10 meV, which can be 
compared with the value of 25-35 meV obtained from fits 
to the spin wave dispersion in the ordered phase. 

With regard to the 2D nature of the spin fluctuations 
above Tjv, we can estimate the relative strength of the 
interlayer coupling in the paramagnetic phase through 
the parameter rjc- The resulting exchange ratio Jd J2 — 
0.1 ± 0.01, is lower than the ratio Jc/J2 ~ 0.15 - 0.25 
in the ordered AFM stateii, implying that the fluctua- 
tions have a more two-dimensional character above Tjv. 
This result is consistent with angle-resolved photoemis- 
sion spectroscopy (ARPES),^? which showed that the 
band structure dispersion of CaFe2As2 is quasi-2D above 
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FIG. 10: Schematic drawing of the anisotropic scattering in the Fe 
layer at Qafm- The inset shows roughly the extent of spin corre- 
lations in the layer which result in anisotropic scattering at Qafm • 
Correlations transverse to Qafm a-re shorter due to the frustrating 
effect of antiferromagnetic Ji interactions on ferromagnetic bonds 
(shown in red), whereas Ji strengthens correlations along the an- 
tiferromagnetic bonds extending in the longitudinal direction. 



T/v but 3D below Tjv. 

At last, we would like to mention that the existence 
of anisotropic in-plane correlation lengths has a sim- 
ple physical interpretation within the frustrated Ji — J2 
model. The magnetically ordered phase consists of di- 
agonal stripes with the magnetic moments ordered fer- 
romagnetically along one diagonal and antiferromagneti- 
cally along the correspondent orthogonal diagonal, as il- 
lustrated in Fig. [TUl The direction longitudinal to Qafm 
points along chains of antiferromagnetic spins while the 
transverse direction points along ferromagnetic chains. 
Since Ji and J2 are antiferromagnetic couplings, the 
nearest-neighbor bonds along the transverse direction are 
frustrated, leading to a shorter correlation length when 
compared to the correlation length associated to the un- 
frustrated bonds along the longitudinal direction. Thus, 
the data are consistent with nanoscale regions of dynamic 
AFM stripe/nematic correlations. 



B. Summary 

In conclusion, paramagnetic excitations in CaFe2As2 
(T5 = T/v = 172 K) have been investigated in the tem- 
perature range from 180 K 1.05rAr) up to 300 K 
(l.ST/v). Above T/v — 172 K, the spin gap collapses 
and antiferromagnetic spin wave scattering is replaced 
by in-plane anisotropic short-ranged AFM correlations 
that extend to high energy [fiuj <60 meV) and tempera- 
ture (r < 300 K). The paramagnetic excitations are ob- 
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served only around the low temperature magnetic zone 
center Qafm- The in-plane anisotropy of the spin corre- 
lation length has a natural explanation in terms of short- 
ranged magnetic interactions in the Ji — J2 model. Thus, 
the paramagnetic fluctuations correspond to the emer- 
gent stripe AFM ordered structure and are suggestive of 
nematic correlations. The magnetic correlations above 
T/v have an anisotropic 2D character and bear a close re- 
semblance to the paramagnetic fluctuations observed in 
the superconducting compositions. 
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